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I. INTRODUCTION 

The current generation of relativistic heavy-ion collision ex- 
periments should exceed the energy density necessary for the 
formation of a quark-gluon plasma. It is therefore necessary to 
have a quantitative theoretical framework which can be used 
to calculate the properties of a quark-gluon plasma. The usual 
line of reasoning is that since QCD is asymptotically free, its 
running coupling constant Oa becomes weaker as the temper- 
ature increases and therefore the behavior of hadronic matter 
at sufficiently high temperature should be calculable using 
perturbative methods. Unfortunately, a straightforward per- 
turbative expansion in powers of Os does not seem to be of any 
quantitative use even at temperatures many orders of magni- 
tude higher than those achievable in heavy-ion collisions. 

The problem can be seen by looking at the perturbative 
expansion of the free energy .7-" of a quark-gluon plasma, 
whose weak-coupling expansion has been calculated com- 
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where \i is the renormalization scale, oia — as{fi) is the run- 
ning coupling constant in the MS scheme, and we have set 
Nc = 3. The coefficient of al log has recently been com- 
puted Q; however, since there are unknown perturbative and 
non-perturbative contributions at 0{al) we do not include 
terms higher than 0{as^^) in Eq. 

In Fig.Q the free energy with Nf = 2 is shown as a function 
of the temperature T/Tc, where Tc is the critical temperature 
for the deconfinement transition. In the plot we have scaled 
the free energy by the free energy of an ideal gas of quarks 
and gluons which for arbitrary Ac and A/ is 
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and ai^^ are shown as bands that correspond to varying the 
renormalization scale, /i, by a factor of two around the cen- 
tral value n = 27rr. As successive terms in the weak-coupling 
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FIG. 1: The perturbative free energy of QCD for iV/ = 2 
massless quarks as a function of T/Tc. The weak-couphng 

3/2 2 5/2 

expansions through orders a^, Os , ctg, and are shown 
as bands that correspond to varying the renormalization scale 
by a factor of two around 2t:T. Also shown is a lattice 
estimate by Karsch et al. 0| for the free energy. The band 
indicates the estimated systematic error of their result which 
is reported as (15±5)%. 



expansion are added, the predictions change wildly and the 
sensitivity to the renormalization scale grows. It is clear that 
a reorganization of the perturbation series is essential if per- 
turbative calculations are to be of any quantitative use at 
temperatures accessible in heavy-ion collisions. 

The free energy can also be calculated nonperturbatively 
using lattice gauge theory 0. The thermodynamic functions 
for pure-glue QCD have been calculated with high precision 
by Boyd et al. 0. There have also been calculations which 
include dynamical quarks 0,13. In Fig. 1 we have included 
the latest lattice estimate of Karsch et al. ■'7j for the free 
energy for Nf = 2 flavors of light quarks. The band indi- 
cates the estimated systematic error of their result which is 
reported as (15±5)%. Note that the quarks in the simulations 
do have non-zero masses and that extrapolation to zero quark 
mass would require significant computing time. Due to the 
difficulty associated with the inclusion of light/massless dy- 
namical quarks on the lattice it is therefore desirable to have 
analytic methods which can be used to estimate the thermo- 
dynamic functions. 

The only rigorous method available for reorganizing per- 
turbation theory in thermal QCD is dimensional reduction to 
an effective 3-dimensional field theory [9lllf]|. The coefficients 
of the terms in the effective lagrangian are calculated using 
perturbation theory, but calculations within the effective field 
theory are carried out nonperturbatively using lattice gauge 
theory. Dimensional reduction has the same limitations as 
ordinary lattice gauge theory: it can be applied only to static 
quantities and only at zero baryon number density. Unlike in 
ordinary lattice gauge theory, light dynamical quarks do not 
require any additional computer power, because they only 
enter through the perturbatively calculated coefficients in the 
effective lagrangian. This method has been applied to the 
Debye screening mass for QCD [lo| as well as the pressure 
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There are some proposals for reorganizing perturbation the- 
ory in QCD that are essentially just mathematical manip- 
ulations of the weak coupling expansion. The methods in- 
clude Pade approximates 11], Borel resummation |l^ . and 
self-similar approximates 13j. These methods are used to 
construct more stable sequences of successive approximations 
that agree with the weak-coupling expansion when expanded 
in powers of Us- These methods can only be applied to quan- 
tities for which several orders in the weak-coupling expansion 
are known, so they are limited in practice to the thermody- 
namic functions. 

One promising approach for reorganizing perturbation the- 
ory in thermal QCD is to use a variational framework. The 
free energy F is expressed as the variational minimum of a 
thermodynamic potential Q.{T,as\m?) that depends on one 
or more variational parameters that we denote collectively by 
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A particularly compelling variational formulation is the <1>- 
derivable approximation, in which the complete propagator 
is used as an infinite set of variational parameters T4j. The 
^-derivable thermodynamic potential Q, is the 2PI effective 
action, the sum of all diagrams that are 2-particle-irreducible 
with respect to the complete propagator [lal . The n-loop $- 
derivable approximations, in which Q. is the the sum of 2PI 
diagrams with up to n loops, form a systematically improv- 
able sequence of variational approximations. Until recently, 
"l?-derivable approximations have proved to be intractable for 
relativistic field theories except for simple cases in which the 
self-energy is momentum-independent. However there has 
been some recent progress in solving the 3-loop "l>-derivable 
approximation for scalar field theories. Braaten and Petitgi- 
rard have developed an analytic method for solving the 3-loop 
^-derivable approximation for the massless <f)*' field theory 
[l^ . Van Hees and Knoll have developed numerical meth- 
ods for solving the 3-loop ^-derivable approximation for the 



massive </> field theory 



They also investigated renormal- 



ization issues associated with the ^-derivable approximation. 
These issues have recently been studied in detail by Blaizot, 
lancu, and Reinosa 17^ . 

The application of the "I>-derivable approximation to QCD 
was first discussed by McLerran and Freedman 119(1 . One 
problem with this approach is that the thermodynamic poten- 
tial Q. is gauge dependent, and so are the resulting thermody- 
namic functions. The gauge dependence is the same order in 
as as the truncation error when evaluated off the stationary 
point and twice the order in as when evaluated at the sta- 
tionary point |20|. However the most serious problem is that 
even the application of 2-loop "I>-derivable approximation to 
gauge theories has proved to be intractable. 

The 2-loop "I>-derivable approximation for QCD has been 
used as the starting point for HTL resummations of the en- 
tropy by Blaizot, lancu and Rebhan |21jj and of the pressure 
by Peshier ^2^. The thermodynamic potential n2-ioop is a 
functional of the complete gluon propagator D;j^(P). How- 
ever, in order to make the problem tractable the authors in 
Refs. I21II and [2^ were forced to make a variational ansatz for 
the exact gluon propagator which they took as the HTL gluon 
propagator in the infrared and free in the ultraviolet with an 
aribitrary momentum scale separating the two momentum re- 
gions. Using this ansatz they were able to calculate the QCD 
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thermodynamic functions; however, a first-principles calcu- 
lation of the corrections to their results for gauge theories 
would require the inclusion of exact vertices as well as exact 
propagators thus making the problem intractable. 

The difficulties in calculating quantities using "l>-derivable 
approximations in gauge theories motivates the use of simpler 
variational approximations. One such strategy that involves a 
single variational parameter m has been called optimized per- 
turbation theory (2^. variational perturbation theory [2^. or 
the linear 5 expansion |2^ . This strategy was applied to the 
thermodynamics of the massless (f)^ field theory by Karsch, 
Patkos and Petreczky under the name screened perturbation 
theory 2G|. The method has also been applied t o sp onta- 
neously broken field theories at finite temperature |27ll . The 
calculations of the thermodynamics of the massless field 
theory using screened perturbation theory have been extended 
to 3 loops 28J. The calculations can be greatly simplified by 
using a double expansion in powers of the coupling constant 
and m/T 0. 

HTL perturbation theory (HTLpt) is an adaptation of this 
strategy to thermal QCD pffll |. The exactly solvable theory 
used as the starting point is one whose propagators are the 
HTL quark and gluon propagators. The variational mass pa- 
rameters mo and mq are identified with the Debye screening 
mass and the induced quark mass. The one- loop free energy 
in HTLpt was calculated for QCD in Ref. (30|| and for QCD 
with massless quarks in Ref. |3lll . At this order, the param- 
eters m_D and could not be determined variationally, so 
their perturbative limits were used. The resulting thermo- 
dynamic functions had errors of order as, but the terms of 
order Qs''^ associated with Debye screening were correct. A 
two-loop calculation is required to reduce the errors to order 

In a previous paper we calculated the thermodynamic 
functions of pure-glue QCD to next-to-leading order in 
HTLpt [3^. In that paper we showed that it was possible 
to renormalize the resulting expressions for the thermody- 
namic potential at next-to-leading order using only vacuum 
and mass counterterms and we also showed that the correc- 
tions to the thermodynamic functions in going from leading- 
order to next-to-leading order were small down to tempera- 
tures on the order of 10 Tc. In this paper we calculate the 
thermodynamic functions of QCD to next-to-leading order in 
HTLpt including the contributions from quark and quark- 
gluon interaction diagrams. 

We begin with a brief summary of HTL perturbation theory 
including quarks in Section [H] In Section IIIII we give the 
expressions for the one-loop and two-loop diagrams for the 
thermodynamic potential. In Section If VI we reduce those 
diagrams to scalar sum-integrals. We are unable to compute 
those sum-integrals exactly, so in Section[V]we evaluate them 
by treating mo and rUq as 0(17) quantities and expand them 
in mo/T and niq/T keeping all terms that contribute up to 
0{g^). The dia grams are combined in Section fVIII to obtain 
the final result for the two-loop thermodynamic potential up 
to 0{g^). In Section [Villi we present our numerical results 
for the free energy of QCD at leading and next-to-leading 
order in HTLpt. In Section ITXI we evaluate the free energy in 
the large Nt limit where exact numerical results have been 
obtained [lllll. 

There are several appendices that contain technical details 
of the calculations. In Appendix EI we give the Feynman 
rules for HTL perturbation theory in Minkowski space to fa- 



cilitate the application of this formalism to signatures of the 
quark-gluon plasma. The most difficult aspect of these calcu- 
lations was the evaluation of the sum-integrals obtained from 
the expansion in mo/T and rUq/T. We give the results for 
these sum- integrals in Appendix^ The evaluation of some 
difficult thermal integrals that were required to obtain the 
sum-integrals is described in Appendix Icl 



II. HTL PERTURBATION THEORY 

The lagrangian density that generates the perturbative ex- 
pansion for QCD can be expressed in the form 

-Cqcd = -]f£r{G^,M'"')+ii>l''D^,■li, 

+Cf,f + £ghost + A£qcd . (9) 

The gauge potential is = with generators of 

the fundamental representation of SU(A'^c) normalized so that 
Trft^ = (5"V2- The field strength tensor is G,,^ = - 
dvAi_i — ig[A^, A,j]. In the quark term there is an implicit 
sum over the Nf quark flavors and = — igA^^ is the 
covariant derivative for the fundamental representation. The 
ghost term £ghost depends on the choice of the gauge-fixing 
term £gf. Two choices for the gauge-fixing term that depend 
on an arbitrary gauge parameter ^ are the general covariant 
gauge and the general Coulomb gauge: 

rgf = --^Tr [{d'^A^.f) covariant , (10) 

= -iTr((V-A)^) Coulomb. (11) 

It is also convenient to introduce various invariants associated 
with the representations of the SU(A^c) gauge group. Denot- 
ing the generators of the adjoint representation as (F")*"^ = 
— i/"*"^ and generators of the fundamental representation as 
we define the following group theory factors: 

[F F ] = / / =ca5 , 

TtFT'' = saS"'', 

5^^ = dA, 

TrTT* = sfS''\ 

S" — dp—SFdA/CF- (12) 

With the standard normalization 

CA = SA = Nc, 

dA = iVc-1, 

CF = (iVf - l)/(2iV,) , 

SF = Nf/2, 

dF = N,Nf. (13) 

The perturbative expansion in powers of g generates ul- 
traviolet divergences. The renormalizability of perturbative 
QCD guarantees that all divergences in physical quantities 
can be removed by renormalization of the coupling constant 
Qfs = (?^/47r. There is no need for wavefunction renormaliza- 
tion, because physical quantities are independent of the nor- 
malization of the field. There is also no need for renormaliza- 
tion of the gauge parameter, because physical quantities are 
independent of the gauge parameter. 
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Hard-thermal-loop perturbation theory (HTLpt) is a reor- 
ganization of the perturbation series for thermal QCD. The 
lagrangian density is written as 

-C = (^QCD + -Chtl) +A£htl. (14) 



The HTL improvement term is 
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where in the first term is the covariant derivative in the 
adjoint representation, in the second term is the covariant 
derivative in the fundamental representation, j/*^ = (l,y) is a 
light-like four- vector, and {■ ■ ■)y represents the average over 
the directions of y. The term HSU has the form of the effective 
lagrangian that would be induced by a rotationally invariant 
ensemble of colored sources with infinitely high momentum. 
The parameter mo can be identified with the Debye screening 
mass and the parameter can be identified as the induced 
finite temperature quark mass. HTLpt is defined by treating 
5 as a formal expansion parameter. 

The HTL perturbation expansion generates ultraviolet di- 
vergences. In QCD perturbation theory, renormalizability 
constrains the ultraviolet divergences to have a form that can 
be cancelled by the counterterm lagrangian A£qcd. We will 
demonstrate that renormalized perturbation theory can be 
implemented by including a counterterm lagrangian A£htl 
among the interaction terms in l|14|l . There is no proof that 
the HTL perturbation expansion is renormalizable, so the 
general structure of the ultraviolet divergences is not known; 
however, it was shown in our previous paper 32j that it was 
possible to renormalize the next-to-leading order HTLpt pre- 
diction for the free energy of pure-glue QCD using only a vac- 
uum counterterm and Debye mass counterterm. Here we show 
that when quarks are included it is also possible to renormal- 
ize the resulting expressions using only vacuum, Debye mass, 
and quark mass counterterms. 

The leading term in the delta expansion of the vacuum 
energy, So, counterterm A£o was deduced in Ref. js^] by cal- 
culating the free energy to leading order in 5. The £o coun- 
terterm A£o must therefore have the form 
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To calculate the free energy to next-to-leading order in 5, we 
need the counterterm Afo to order 5 and the counterterms 
Am|, and Am^ to order 5. We will show that there is a non- 
trivial cancellation of the ultraviolet divergences if the mass 
counterterms have the form 



Am£3 = 
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Physical observables are calculated in HTLpt by expanding 
them in powers of 5, truncating at some specified order, and 
then setting 8 = 1. This defines a reorganization of the per- 
turbation series in which the effects of the mf, and terms 



in (list are included to all orders but then systematically sub- 
tracted out at higher orders in perturbation theory by the 
Smjj and Sm^ terms in 1151 . If we set 5 = 1, the lagrangian 
II14II reduces to the QCD lagrangian @. If the expansion in 5 
could be calculated to all orders, all dependence on mo and 
TTiq should disappear when we set 5 = 1. However, any trun- 
cation of the expansion in 5 produces results that depend on 
mo and lUq. Some prescription is required to determine mo 
and ruq as a function of T and Qs. We choose to treat both 
as variational parameters that should be determined by mini- 
mizing the free energy. If we denote the free energy truncated 
at some order in 5 by 0.{T, a^, mD,rnq, 5), our prescription is 



d 
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dm„ 



Q{T,as,mD,mq,5 = 1) = 0, 
Q,{T,as,mD,mq,5 = 1) = 0. 



(19) 



Since Q{T,as,mD,mq,S = 1) is a function of the variational 
parameters mo and ruq, we will refer to it as the thermody- 
namic potential. We will refer to the variational equations 
I19II as the gap equations. The free energy is obtained by 
evaluating the thermodynamic potential at the solution to 
the gap equations 1191 1. Other thermodynamic functions can 
then be obtained by taking appropriate derivatives of with 
respect to T. 



III. DIAGRAMS FOR THE THERMODYNAMIC 
POTENTIAL 

The thermodynamic potential at leading order in HTL per- 
turbation theory for an SU{Nc) gauge theory with Nf mass- 
less quarks is 



i^LO = dA^g + dp^q + AqSo 



(20) 



where J^g is the contribution from each of the color states of 
the gluon: 

^9 = -V {(d-l)loghAT(P)] + logAi(P)} .(21) 

The transverse and longitudinal HTL propagators At (-P) and 
Al{P) are given in 1A.5H and 1A.52II . The quark contribu- 
tion is 



T{p} 



logdet[f- E(P)] 



(22) 



where E(P) is the HTL fermion self-energy. The leading order 
counterterm AqSo was determined in Ref. 1301 
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(23) 



The thermodynamic potential at next-to-leading order in 
HTL perturbation theory can be written as 



f^NLC 



OlO + dA [J'Sg + J~Ag + T gh + J~gci\ 
-\-dASF + T4,qg\ + dfTqct + AiiSq 
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where AiEq and Aim|j are the terms of order S in the vacuum 
energy density and mass counterterms. The contributions 
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FIG. 2: Diagrams contributing through NLO in HTLpt. 
Shaded circles indicate dressed HTL propagators and vertices. 



from the two-loop diagrams with the three-gluon and four- 
gluon vertices are 



xr"""(P,Q,ii)A'"'(P)A^"(Q)A''"(i?) , (25) 
r^"'''"(P,-P,Q,-Q) 

PQ 

xA^''(P)A^"(Q) , (26) 
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where -Q- P. 

The contribution from the ghost diagram depends on 
the choice of gauge. The expressions in the covariant and 
Coulomb gauges are 

= T^'^^^^i^^''^'^'"'^-^) covariant, (27) 

= ^/if W{Q^-Q-nn'^){R'' -R-nn'') 
2 TpQ r 

xA'"'{P) Coulomb . (28) 

The contribution from the HTL gluon counterterm diagram 
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The contributions from the two-loop diagrams with the 
quark-gluon three and four vertices are given by 

•^3,9 = i/^f Tr[r^(P,Q,P)S(Q) 
^ JP{Q} 

xV'{P,Q,R)S{R)]A'"'{P) (30) 

^4,9 = iff'^f Tr[r^''(P,-PQ,0)S(0)] 



xA^''(P) , 

where Tr implies taking the trace over 7-matrices. 
tribution from the HTL quark counterterm is 
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(31) 
The con- 
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Provided that HTL perturbation theory is renormalizable, 
the ultraviolet divergences at any order in S can be cancelled 
by renormalizations of the vacuum energy density £0 , the HTL 
mass parameters mf, and nig, and the coupling constant as. 
Renormalization of the coupling constant does not enter un- 
til order S^. We will calculate the thermodynamic potential 
as a double expansion in powers of mo/T and mq/T, ad g 
including all terms through 5*** order. The Sos term in A£o 
does not contribute until 6**^ order in this expansion, so the 
term of order 5 in A£o can be obtained simply by expanding 
Eq. 12311 to first order in S: 
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(33) 



The remaining ultraviolet divergences must be removed by 
renormalization of the mass parameters ran and m,. We will 
show below that, at order 8, all remaining divergences can be 
removed by the quark and Debye mass counterterms. This 
provides nontrivial evidence for the renormalizability of HTL 
perturbation theory at this order in 5. 

The sum of the 3-gluon, 4-gluon, and ghost contributions is 
gauge invariant. By using the Ward identities, one can easily 
show that the sum of these three diagrams is independent 
of the gauge parameter ^. With more effort, one can show 
the equivalence of the covariant gauge expression with ^ = 
and the Coulomb gauge expression with ^ = 0. In a similar 
manner, it can be shown that the sum of l|3()|l and H31|l is 
independent of ^ within the class of covariant and Coloumb 
gauges, as well as the equivalence of the two with ^ = 0. 



IV. REDUCTION TO SCALAR 
SUM-INTEGRALS 

The first step in calculating the quark contribution to the 
free energy is to reduce the sum of the diagrams to scalar- 
sum-integrals. The leading-order quark contribution can be 
rewritten as 



2^ logP^-2^ log 
T{P} T{P} 



Al-Al 



p2 



where 



Ao{P) = iPo - T^Tp , 
iPo 

As{P) = \p\ + -^[l-Tp] . 

IpI 

The HTL quark counterterm can be rewritten as 
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(34) 

(35) 
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(37) 



6 



We proceed to simplify the sum of USUI and 1311 in Landau gauge. Using the Ward identities IIA.45t and 1A.48II the sum of 
l|3()|l and becomes 

■^3,a+4,« = i/^^ ^ |Ax(P)Tr [r~5(Q)] - AT(P)Tr [r^S{Q)r^S{R')] + Ax(P)Tr [r05(Q)r»5(E')] | , (38) 

where S is the quark propagator, At is the transverse gluon propagator, Ax is a combinaiton of the transverse and longitudinal 
gluon propagators defined in 1A.27II . and R' = Q — P. 
Performing the traces of 7-matrices gives 



A-rAs{Q)As{R) - Ao{Q)AoiR) 



Al{R)~Al{R) . -V ^\ (p.y)2_(Q.y)2 (Q.y) 



+ 



8m^AT(P) / (AoiQ) - As(Q)q-y)(Ao(P) ^ As{R)i-y) 
%{R)-Al{R)\ {Q-Y){R-Y) 



4mlAx{P) / 2AoiR)As{Q)ci-y-Ao{Q)Ao{R)-As{Q)As{R)ci-r \\ 2 4. 

+ AUR)-AUR) \ (Q-Y){R-Y) /_ f + '"'^ ' ^^'^ 



where Ao and As are defined in llA.40t and llA.41l i. respectively. 

r 



V. HIGH-TEMPERATURE EXPANSION 

The free energy has been reduced to scalar sum-integrals. 
If we tried to evaluate the 2-loop HTL free energy exactly, 
there are terms that could at best be reduced to 5-dimensional 
integrals which would have to be evaluated numerically. We 
will therefore evaluate the sum-integrals approximately by ex- 
panding them in powers of mn/T and niq/T. We will carry 
out the expansion to high enough order to include all terms 
through order if mo and are taken to be of order g. 

The free energy can be divided into contributions from hard 
and soft momenta. We proceed to calculate the hard-hard 
and hard-soft contributions. There is no soft-soft contribution 
since one of the momenta in the loop is always fermionic and 
therefore hard. 



1287r^ 
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The hard contribution from the HTL counterterm I29II was 
given in p3| and reads 



27r^ 

- - 7 + 27 + — 
e 3 



(41) 



The sum-integrals over P involve two momentum scales 
and T. Since Pq — {2n+l)nT, the momentum is always hard. 
We can therefore expand in powers of rrig. To second order 
in nig, we obtain 



A. One-loop sum-integrals 

The one-loop sum-integrals include the leading quark con- 
tribution (I22II and the HTL quark counterterm 1321 1. The 
leading order free energy must be expanded to order to 
include all terms through order if is taken to be of order 

g- 



1. Hard contributions 

The hard contribution from the LO gluon term 12111 was 
given in |33 and reads 



' 45 24 
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log - 4mq 
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1 

p2 



p2 • p2p2'^P p2p2 'JP) 



(42) 



Note that the function Tp cancels from the term. The 
values of the sum-integrals are given in Appendix B. Inserting 
those expressions, the hard quark contributions to the free 
energy reduce to 
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(43) 



Note that this contribution is finite and so the leading order 
counterterm AqEo is the same as in the pure-glue case. The 



HTL quark counterterm is given in 13711 . Expanding this term 
to second order in yields 




The values of the sum-integrals are given in Appendix B. In- 
serting those expressions, the hard contributions to the HTL 
quark counterterm reduce to 

^;^^=-i-^^^-^(-^-6)m^ (45) 

Note that the first term in Eq. 14511 cancels the order-e° term 
in coefficient of in Eq. 1431 



2. Soft contributions 



There are no soft contribution from the leading-order quark 
term l|84|l or from the HTL quark counterterm H37|l . 



B. Two-loop sum-integrals 

The sum of the two- loop sum-integrals is given in l).'-{8^ . 
Since these integrals have an explicit factor of g^, we need only 
to expand the sum-integrals to order nigmD /T'^ and m^/T'' 
to include all terms through order g^. 

The sum-integrals involve two momentum scales , m d 
and T. In order to expand them in powers of.., we separate 
them into contributions from hard loop momenta and soft 
loop momenta. This gives two separate regions which we will 
denote (hh) and {hs). In the (hh) region, all three momenta 
P, Q, and R are hard. In the (hs) region, two of the three 
momenta are hard and the other soft. 



The soft contribution comes from the Po = term in the 
sum-integral. At soft momentum P = (0,p), the HTL self- 
energy functions reduce to Ht(-P) ~ and I1l{P) ~ rnjj. 
The transverse term vanishes in dimensional regularization 
because there is no momentum scale in the integral over p. 
Thus the soft contribution comes from the longitudinal term 
only. 

The soft contribution to the leading order free energy 12111 
was given in Ref. [s^l and reads 



Contributions from the (hh) region 



For hard momenta, the self-energies are suppressed by 
m\,/T'^ or niq/T^ relative to the propagators, so we can ex- 
pand in powers of Ht, Hl, and E. 

The (hh) contribution from 125I I- 127II was given in Ref. [33| 
and reads 



1 

12^ 



1 



8 



mlT. 



(46) 



The soft contribution to the HTL gluon counterterm l|29|l 
given in Ref. j32] and reads 



(hh) 



ag + ig + gh 



12 3tt 



y4 



7_ 

96 



mlT\ 



4-4.621 



3tv 



(48) 



r, 



get 



(47) 



The (hh) contribution from l|3()|l and Ij81|l can be written as 



'{PQ} 



P{Q} 
Ad 



P^Q^ 



+ 2m 



2 2vr 

og 21 



'P{Q} 
2d PQ 



^ 2 2^ [S-d 1 



'{PQ} 
+2my{d~l) 



+2my{d-l) 



{PQ) 

d+1 1 

d-1 P^Q^r^ d - 1 P2Q2^4 _ 1 P2Q2r4 

2d PQ d + 2 1 



p^P'^Q'^ 
Ad 



d- 2 



(P2)2Q2 d-lp2p2Q2 



+ ■ 



)2P2 d_lp2Q2^4 d-lp2Q2^2 d - 1 p2g2r4 d-1 P^Q^r^R'^ 



r 1 



4- 



p2_r2 



PQ} 



{PQ) 



[P'^QlQ^ P^^QlR^ 
d + 3 1 2 



Tq +2m,ff (d- 1 



4 

JP{Q} 



P2(Q2)2 P2Q2Q: 



+ 



d-lP2Q2i^2 P2(Q2)2 q2p2Q2Jl2 



(49) 



Using the expressions for the sum-integrals in Appendix B, this reduces to 
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2. The (hs) contribution 



In the (hs) region, the momentum P is soft. The momenta Q and R axe always hard. The function that multiphes the soft 
propagator At(0, p) or Ax(0,p) can be expanded in powers of the soft momentum p. In the case of At(0, p), the resulting 
integrals over p have no scale and they vanish in dimensional regularization. The integration measure scales like m^, the 

soft propagator Ax(0, p) scales like l/m|,, and every power of p in the numerator scales like m_D. 
The (hs) contribution from 12511 - 1271 1 was given in Ref. 13211 and reads 



.(hs) 



3ff+434 



11 
32 



(51) 



The only terms that contribute through order g^m'^T and m^mog'^T from 13UII and 1311 1 are 



(hs) 



3qg+4qg 



g^T 



1 



p-' + mj, 









m 


4 


3 







(3 + d) jT^^^ttt + -- 



(52) 



In the terms that are already of order g^m^T, we can set 
R = —Q. In the terms of order g^moT^, we must expand 
the sum-integrand to second order in p. After averaging over 
angles of p, the linear terms in p vanish and quadratic terms 
of the form p'p^ axe replaced by p'^5^-' /d. We can set p^ — 
—m\,, because any factor proportional to p^ + mjj will cancel 
the denominator of the integral over p, leaving an integral 
with no scale. This gives 



(hs) 



3<js+4gg 



t ^'i Cts 

2e 



+ l + 27 + 41og2 



(53) 



3. The (ss) contributions 



The (ss) contribution from 125II - 127I I was given in Ref. |32 
and reads 



{ss) 



39+49+9/1 



_3_ 

16 



+ 3 



37r 



There is no (ss) contribution from the diagrams involving 
fermions. 



VI. HTL-IMPROVED THERMODYNAMICS 



The free energy at second order in HTL perturbation theory 
defines a function Q{T,as,mD,rnq,S = 1). We will refer to 
this function as the thermodynamic potential. To obtain the 
free energy F{T) as a function of the temperature, we need 
to specify a prescription for the mass parameter mo as a 
function of T and as . 



S expansion and then to next-to-leading order. 



A. Leading order 

The complete expression for the leading order thermody- 
namic potential is the sum of the contributions from 1-loop 
diagrams and the leading term 1231 in the vacuum energy 
counterterm. The contributions from the 1-loop diagrams, 
including all terms through order g^, is the sum of 14UL I43L 
and lESJ 



loop 



= —dA 



+30ml + ^ 



n I 
45 
1 
e 

dp 



7 dp 15 , 2 nnd-F . 



5,2 
<3 



21og|-7 + 27-f — 



. 4 
mo 



-60^ (.^_ 6)^^} 



(55) 



where rho, rhq and jl are dimensionless variables: 

mo 
= 2^' 
m„ — — — , 

- = M 
^ 2nT ' 

Adding the counterterm 12311 . we obtain the thermodynamic 
potential at leading order in the delta expansion: 



(56) 
(57) 
(58) 



= -dA 



45 
45 
4 

dp 



r 7 dF 15 ^ 2 ^r^dp 



- 2 

rUq 



^"§2-2+^+^ 



. 4 



-60^ (.^-6)^^} 



(59) 



VII. THERMODYNAMIC POTENTIAL 



B. Next-to-leading order 



In this section, we calculate the thermodynamic potential 
0,{T, as,mD,mq, 5 = 1) explicitly, first to leading order in the 



The complete expression for the next-to-leading order cor- 
rection to the thermodynamic potential is the sum of the 
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contributions from all 2-loop diagrams, the quark and gluon 
counterterms, and renormalization counterterms. The contri- 
butions from the 2-loop diagrams, including all terms though 



order is the sum of (gSJ, (HOJ, IHJ, (HSJi and multi- 
plied by the appropriate group structure constants listed in 



fig 



2 — loop 



, Tv^r Q, 5 / 5 \ , „ 55 

= ~dA—z < - T C-4 + -SF + 15(CA + SFjmo — 

45 TT I 4 V 2 / 8 



(cA - Yl^F^ +4 log 



'72 



^ - 0.4714 SF 


. 2 45 r 


mo —sf 







i 4 log I +8.96751 



+ 180 SFrhorh^ 



+ - 



165 



CA~^SF^ +41og^ -21ogmD) +CA +27) - ^SF(l + 27 + 41og2) 



ml, 



The HTL gluon counterterm is the sum of 14H and 1471 



, TT^r 15 „ 2 ... 3 45 / 1 „, A ^ „ 27r , 4 
rfA^g- < -jrhl - 45m?, _ _ f - 2 log | - 7 + 27 -f — ) ml. 



The HTL quark counterterm is given by 11451 



fiqct = -dF^^{30m^-H20(7r^-6)m*} 



(60) 



(61) 



(62) 



The ultraviolet divergences that remain after these 3 terms are added can be removed by renormalization of the vacuum energy 
density £0 and the HTL mass parameter mo- The renormalization contributions at first order in 5 are 



2d 2 ^ 

Af2 = Aifo + Aimj,- — ^-f^LO + Aim,- — :rQ.i,o 
ami Hin^ 



' dmi 



(63) 



Using the results listed in Eqs. 11171 . 11181 . and 11331 the complete contribution from the counterterm at first order in 5 is 



ASl = -dA 
165 



TT^r* I 45 



45 Il^^'^ + V 



f i + 21ogf+2i^+2).^ 



(ca-^sf) (i + 21og ^-21ogmi,+ 2) mi, + ySF (^i + 2 + 2 log | - 2 log 2 -f 2^^^ m^ | 



(64) 



Adding the contributions from the two- loop diagrams in 1601 . the HTL gluon and quark counterterms in 1611 and 1621 . the 
contribution from vacuum and mass renormalizations in 1641 , and the leading order thermodynamic potential in 1591 we obtain 
the complete expression for the QCD thermodynamic potential at next-to-leading order in HTLpt: 

fiNLO = -'^^^ l+4d7-15"^o-T (^°S2 -2+7+y )m^+60— - 6) m. 



TT 



001) - i-SF (log I -2.337) \ml 



-45 Sf 



5 / 5 \ , ^ , , , 55 f /, u 36 , . 

4 \^ 2^^) ^ ^'^^ + SFjmo - —{ CA (log-^ - Yi - 2 

(log I + 2.192) m^ + 1^ jcA (log I + A + ^) _ ±SF (log I - i + 7 + 2 log 2)|mi, + 180 SFmom\ | 



.(65) 



C. Gap Equation 



Onlo with respect to each parameter taken holding the other 
constant vanishes 



The quark and gluon mass parameters, mq and mo, are 
determined variationally by requiring that the derivative of 



d 

T. ilNLo(r, as, ?rii3, m,) — 

om„ 



(66) 
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10"' 10 
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FIG. 3: Numerical solution of gap equations for (a) tud I69II 
and (b) 1681 1 as a function of as(27rr) for Nc = 3 and 
Nf = 3. The shaded band corresponds to varying the renor- 
malization scale p by a factor of two around ^ — 2-kT. 



f^NLO (T, Qs , mo , rriq 



(67) 



The first equation above results in the following equation for 



dA 



{n 



6)mq = 



3 I log ^+2.192) -12mD 



(68) 



In the limit of small Os the above gap equation does not 
go to the perturbative limit for the quark mass which is 
"iq.pert = cf as/Svr. The fact that does not go to the 
perturbative value in the small limit is due to the fact 
that the perturbative limit of the quark gap equation results 
from terms which are 0(aa) and these terms are not included 
completely at NLO in HTLpt. One might hope that going to 
the next order in HTLpt would cure this problem; however, 
this will in fact not happen since the fermion sector is infrared 
safe and therefore only even powers of rhq will appear at each 
order. At NNLO all terms contributing at 0{ag) at NLO will 
be replaced by explicit powers of as and all rhq dependence 
will be pushed up to 0(as). This behavior will persist at all 



orders in HTLpt so that at any order the weak-coupling limit 
of the gap equation quark mass will be scale dependent. In or- 
der to circumvent this problem we can consider other possible 
prescriptions for choosing rhq which include requiring that rhq 
be equal to its perturbative value for all as or requiring that 
rhq be proportional to rho with the proportionality constant 
fixed in the weak-coupling limit. 

Performing the derivative with respect to m_D while holding 
niq fixed results in the following gap equation for mo 



45m|) 



1 + 



l°S2-2+'^+T 



a I ^ , , 55 



CA 



(log 



A 
2 



36 
11 



logmu — 3.637 



11 



11 



sf 



Sf 



log 



lo. 



2 



2.337 

1 
2 



rho + ■ 



- 7 + 2 log 2 



495 



CA (log ^ + 



22 



+ 7 



m|) + 180 SFihg 



(69) 



The last term in Eq. 16911 proportional to rhq can be written 
in terms of rhn using Eq. 1681 . In Fig. |3]we plot the solutions 
to the gap equations for mo and rUq for Nc = 3 and Nf = 
2. The solution for mo goes to the perturbative value in 
the limit of small as, decreases below the perturbative value 
as as increases, and becomes larger than the perturbative 
value at as ~ 0.11. The solution for rriq does not go to the 
perturbative value in the limit of small as and is instead scale 
dependent even at lowest order as dicussed above. As as 
increases the value of remains very flat regardless of the 
scale, changing significantly only near as ~ 0.10. 



VIII. FREE ENERGY 

The free energy is obtained by evaluating the leading and 
next-to-leading order thermodynamic potentials, 1591 and 
16511 ■ at the solution to the gap equations 1681 1 and 16911 . In 
Fig. 1^ we plot the leading and next-to- leading order HTLpt 
predictions for the free energy of QCD with A'^c = 3 and 
Nf = 2. We have studied the alternative prescriptions for 
the quark mass discussed in the previous section and find 
that the NLO free energy obtained using these prescriptions 
is numerically indistinguishable from that obtained using the 
quark gap equations. As can be seen from this figure the cor- 
rections in going from LO to NLO are small over the entire 
temperature range, especially when compared to convergence 
of the perturbative result. Additionally, the scale variation 
of the NLO HTLpt result for the free energy is much smaller 
than the LO showing that the partial resummation of the 
scale dependent logarithms reduces the scale variation of the 
final results significantly. 

However, as was the case in pure-glue QCD f33|, the results 
seem to lie significantly above the lattice data which is avail- 
able below 5Tc. There are several reasons for why HTLpt 
might fail to describe the lattice data in this temperature 
range. One is that the hard modes are not resummed prop- 
erly within HTLpt and that a description using a "I'-derivable 
approach which explicitly separates the hard and soft modes 
as done in Ref. 21] is better. A second is that HTLpt dis- 
cards some important physics like topological modes or the 
Zn symmetry of QCD near the phase transition. 
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1000 
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FIG. 4: LO and NLO HTLpt predictions for the free energy of 
QCD with Nc — 3 and Nf = 2 together with the perturbative 
prediction accurate to g''. The shaded bands correspond to 
varying the renormaUzation scale ^ by a factor of two around 
= 2itT. Also shown is a lattice estimate by Karsch et 
al. i for the free energy. The band indicates the estimated 
systematic error of their result which is reported as (15±5)%. 



A third possibility is that the expansion in mo/T and 
niq/T breaks down at these temperatures. Numerically, 
mn/T ~ 1.2 and rUg/T ~ 0.5 at 5T^ and mo/T ~ 1.6 and 
vriq/T ~ 0.6 at 2Tc, which casts doubt on the applicability 
of the expansion in this temperature range. However, in the 
case of pure-glue we have been able to compare the LO HTLpt 
result expanded to 0(m|j) with the non-truncated LO expres- 
sion which is accurate to all orders in m_D and find that the 
expansions converge very rapidly. Numerically we find that at 
niD/T = 5 truncations of the LO order result accurate to m|j 
and rh% reproduce the exact result to 5% and 0.2%, respec- 
tively. There have also been studies of the convergence of the 
mass expansions of the three-loop free energy for a massless 
scalar field theory using screened perturbation theory 29] and 
the ^-derivable approach fl^ which demonstrated that mass 
expansions also converge very rapidly at NLO and NNLO. 
This gives us some confidence that the truncated NLO solu- 
tions are numerically reliable. 
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FIG. 5: NLO HTLpt prediction for the O(iV^) contribu- 
tion to the free energy, the numerical result of Ref. jsji 
and the perturbative prediction accurate to as a func- 
tion of 3cff(MDR) = \/s75(mdr) = 2tv ^/7faJJjtYm) where 
Mdr = TTe~''T. Dots indicate the point at which there is 
no longer a real- valued solution to the gap equation for mo. 
In (a) the renormalization scale n is varied by a factor of e 
around /idr. In (b) the renormalization scale /i is varied by a 
factor of 2 around 2nT. In both (a) and (b) the perturbative 
g^ result is evaluated at the central scale. 



IX. LARGE Nf 

In the limit that Nf is taken large while holding g'^Nf fixed 
it is possible to solve for the O(N^) contribution to the free 
energy exactly HHil. In Fig.|Slwe plot the NLO HTLpt pre- 
diction for the 0{Nf) contribution to the free energy along 
with the numerical result of Ref. and the perturbative 
prediction which is accurate to 0{as^^). In Fig. |S|we plot 
the NLO HTLpt prediction for mo in the large Nf limit and 
the exact numerical result [ssll . The HTLpt predictions for 
both the free energy and the Debye mass seem diverge from 
the exact result around g^s ~ 2 regardless of the scale which 
is chosen; however, for both quantities, choosing the scale to 
he fi = pdr = ■ne~~'T seems to provide a reasonable repro- 
duction of the exact results. This result is comparable to the 
performance of the ^-derivable approach in the large Nf limit 




X. CONCLUSIONS 

In this paper we have extended our previous HTLpt cal- 
culation of the thermodynamic functions in pure-glue QCD 
to include the contribution of Nf massless quarks. We have 
presented results for the leading- and next-to-leading-order 
HTLpt predictions for the QCD free energy for arbitrary Nf . 
Using the NLO HTLpt expression for the thermodynamic po- 
tential we were able to find variational solutions for both the 
quark and gluon mass parameters allowing a first-principles 
prediction of the QCD free energy. As in the case of pure- 
glue we find that the NLO HTLpt prediction lies significantly 
above the available lattice data below 5 Tc\ however, the prob- 
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Exact Result 
NLO HTLpt 




NLO HTLpt agrees with the exact result only out to ^cff ~ 2 
and has large scale dependence after this point. The large 
scale dependence is not surprising given the fact that in the 
large Nf limit the running of the coupling constant is driven 
by the presence of the Landau singularity and even the exact 
results are sensitive to this beyond gcff ~ 5. The poor per- 
formance of NLO HTLpt, however, is comparable to recent 
large Nf predictions within the ^-derivable approach. The 
failure of both approaches to agree better with the exact re- 
sult for large values of g^g is an indication that a description 
of strongly-coupled QCD thermodynamics solely in terms of 
HTL quasiparticles is perhaps inappropriate. However, it is 
possible that the physics of large- A'^/ QCD is so different from 
that of QCD with a small number of flavors that it cannot 
serve as a definitive testing ground for the a ppli cability of the 
quasiparticle approach to the physical case |3q|. 
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NLO HTLpt 



fI = 7lT 

|.L = 2tiT 



|i = 47iT """----,„ 
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APPENDIX A: HTL FEYNMAN RULES 
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1 2 3 4 5 

geff(^DR) 

FIG. 6: NLO HTLpt prediction for mo in the large Nf 
limit and the exact numerical result [ssll as a function of 

ffeff(AtDR) = ^/sJg{^iDR) = 2tv y/ Sfasi^Bn) where /^dr = 
Tve 'T. In (a) the renormalization scale jj, is varied by a 
factor of e around fion- In (b) the renormalization scale /i is 
varied by a factor of 2 around 2nT. 



In this appendix, we present Feynman rules for HTL per- 
turbation theory in QCD. We give explicit expressions for the 
propagators and for the quark-gluon 3 and 4 vertices. The 
Feynman rules are given in Minkowski space to facilitate ap- 
plications to real-time processes. A Minkowski momentum is 
denoted p — {po, p), and the inner product is p-q — poqo~p-q. 
The vector that specifies the thermal rest frame is n = (1, 0). 



1. Gluon Self-energy 



lem of oscillation of successive approximations and large scale 
dependence of the perturbative results is eliminated by using 
this reorganization. 

The failure of HTLpt to describe the lattice data in this 
region could be attributed to the failure of the expansion per- 
formed in mo and rhq; however, a study of the convergence 
of the truncated LO expressions to a numerical evaluation of 
the exact LO expression shows that these expansions converge 
very rapidly. Therefore, we are steered towards the conclu- 
sion that a systematic description of QCD thermodynamics 
using HTLpt is not appropriate below 5Tc. The "l>-derivable 
approach seems to agree better with the lattice data in this 
range so perhaps HTLpt is not resumming the hard modes 
properly and an explicit separation of hard and soft scales is 
required. However, we should point out that some authors 
believe that a description of QCD thermodynamics near the 
phase transition in terms of Polyakov loops is necessary |37|. 

We have also compared the NLO HTLpt free energy and 
Debye mass with exact results which are available in the large 
Nf limit. This comparison shows that, in the large Nf limit. 



The HTL gluon self-energy tensor for a gluon of momentum 
p is 

n^-{p) = ml [T'"'(p, -p) - . (A.l) 

The tensor T'"'(p, q), which is defined only for momenta that 
satisfy p + q = 0, is 

T^-'{P,-P) = (y^y''^) . (A.2) 

The angular brackets indicate averaging over the spatial di- 
rections of the light-like vector y = (l,y)- The tensor T*^" is 
symmetric in fx and f and satisfies the "Ward identity" 

p^T^'^ip, -p) = p-n n . (A. 3) 

The self-energy tensor H^" is therefore also symmetric in fi 
and u and satisfies 

p^n'^-'ip) = 0, (A.4) 
g^.n'^^p) = -ml. (A.5) 
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The gluon self-energy tensor can be expressed in terms 
of two scalar functions, the transverse and longitudinal self- 
energies Ut and IIi,, defined by 

nrip) = j^{5'' -pp')n^'ip), (A.6) 



(A.7) 



where p is the unit vector in the direction of p. In terms of 
these functions, the self-energy tensor is 



The four-vector rip is 



Tip rip 



n-p 



(A.8) 

(A.9) 
(A.IO) 

(A.ll) 



and satisfies p-Up = and rip = 1 — (n-p)^ /p^ . The equa- 
tion (IA.5I I reduces to the identity 



(d-i)nT(p) + — nL(p) = ml 



(A.12) 



We can express both self-energy functions in terms of the 
function T°° defined by (TOl : 



Tlrip) 



{d~l)nl 



T"«(p,-p)-l + n^] , (A. 13) 



(A.14) 



In the tensor T^'^{p,—p) defined in ljA.2fl . the angular 
brackets indicate the angular average over the unit vector y. 
In almost all previous work, the angular average in llOl has 
been taken in d = 3 dimensions. For consistency of higher 
order radiative corrections, it is essential to take the angular 
average in d = 3 — 2e dimensions and analytically continue to 
d = 3 only after all poles in e have been cancelled. Expressing 
the angular average as an integral over the cosine of an angle, 
the expression for the 00 component of the tensor is 



w{<l) . 2^-£ PO 

' dc (1 — c ) 



where the weight function w(e) is 

r(2-26)^2, _ r(| 



wie) 



v-^ii-e) r(|)r(i-e) 



,(A.15) 



(A.16) 



The integral in l)A.15|l must be defined so that it is analytic 
at po = 00. It then has a branch cut running from po — — |p| 
to po = -|-|p|- If we take the limit e ^ 0, it reduces to 

^00/ N Po , po + IpI /a 1'7^ 

T (P,-P) = ^^og — , (A.17) 

2|p| po - IpI 

which is the expression that appears in the usual HTL self- 
energy functions. 



The Feynman rule for the gluon propagator is 

ir'A^^p) , (A.18) 

where the gluon propagator tensor A^^ depends on the choice 
of gauge fixing. We consider two possibilities that introduce 
an arbitrary gauge parameter ^: general covariant gauge and 
general Coulomb gauge. In both cases, the inverse propagator 
reduces in the limit ^ ^ cxa to 

A-^(p)'"' = 



2 LIV , LI U 

'P g +p'^p 



W'^ip) . (A. 19) 



This can also be written 

At(p) nf,/\L(p) 

where At and Al are the transverse and longitudinal prop- 
agators: 

1 



At(p) 
Al(p) 



p2-nT(p) ' 
1 



-n2p2 + ni(p) ■ 

The inverse propagator for general ^ is 

A~^(p)'"' = A^^(p)'"' — -^p'^p'^ covariant 



(A.21) 
(A.22) 

(A.23) 



= A-i(p) 

The propagators obtained by inverting the tensors in IIA.241 
and 1IA.23I1 are 



— - (p*" — p-n n^) (p" — p-n n") 

Coulomb . (A. 24) 



A'^-(p) 



~ATip)Tp^'' + ALipKn; - ^ 



(p2)2 

covariant , (A. 25) 



-AT(p)r;'^ + Ai(p)nV-C ^^^"^ 



(ripP^) 
Coulomb . (A. 26) 

It is convenient to define the following combination of prop- 
agators: 



Ax{p) = Ai(p) + — At(p) 



(A.27) 



Using llX?T2l . iTOn . and 11X221 . it can be expressed in the 
alternative form 

Ax(p) = [ml, -dnT(p)] Ai(p)AT(p) , (A.28) 

which shows that it vanishes in the limit mo — > 0. In the 
covariant gauge, the propagator tensor can be written 

A''-'(p) = [-At(p)(7^'' + Ax{p)n^n] 

-^^Aa-(p) (p^n'^ + n^p") 



p^ 



+ 



At(p) + 



{n-pf 



Ax(p)-4 



(A.29) 



This decomposition of the propagator into three terms has 
proved to be particularly convenient for explicit calculations. 
For example, the first term satisfies the identity 

[-At(p)3m-' + Ax{p)n^n^] A'^ipY^ = 

A PmP^ , n-p Axjp) X 
„5 + „2_2 A_/„^^'^'"p■ 



p^ 



n|p2 Al(p) 



(A.30) 
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2. Quark Self-energy 

The HTL self-energy of a quark with momentum p is given 



by 



where 



p-y 



(A.31) 



(A.32) 



Expressing the angular average as an integral over the cosine 
of an angle, the expression is 



■w{e) 



2 /_i ■ ' Po-|P|c 



(A.33) 



The integral in IIA.33|l must be defined so that it is analytic 
at po = oo. It then has a branch cut running from po — — |p 
to po ~ +\p\- In three dimensions, this reduces to 



S(P) = ^70 log ^ 
2|p| Po - IpI 



+ t47 ■ P 1 - 7^ log ^ . (A. 34) 

P V 2 p PO - P ' 



3. Quark Propagator 



The Feynman rule for the quark propagator is 



iS^^Sip) 



The quark propagator can be written as 

1 



(A.35) 



S{P) = 



i> ~ s(p) 



(A.36) 



where the quark self-energy is given by IIA.3H . The 
quark propagator can be written as 

S~\p)^i,-^{p) . 

This can be written as 

S~\p) =4(p) , 

where we have organized ^a(p) and As(p) into: 

A^{p) = (Ao(p),As(p)p) . 

The functions Ao(p) and As{p) are defined 



^o(fi) = Po , 



""9/7- 

Pi) 



As{p) = lpi + ^[l-Tp] 



(A.37) 
(A.38) 
(A.39) 

(A.40) 
(A.41) 



4. Quark-gluon vertex 

The quark-gluon vertex with outgoing gluon momentum p, 
incoming fermion momentum q, and outgoing quark momen- 
tum r, Lorentz index /j, and color index a is 



Ta (P, <l,r) ^ gta (7'' - m^f^ (p, g, r )) 



(A.42) 



The tensor in the HTL correction term is only defined for 
p — q + r = Q: 



f''{p,q,r) = (y' 



i 



q-y r-y 



(A.43) 



This tensor is even under the permutation of q and r. It 
satisfies the "Ward identity" 

p^f^{p,q,r)^t^{q)~t^{r) . (A.44) 

The quark-gluon vertex therefore satisfies the Ward identity 

p^r^ip,q,r) = S-\q)^S-'ir) . (A.45) 

5. Quark-gluon four-vertex 

We define the quark-gluon four-point vertex with outgo- 
ing gluon momenta p and q, incoming fermion momentum r, 
and outgoing fermion momentum s. Generally this vertex has 
both adjoint and fundamental indices, however, for this cal- 
culation we will only need the quark-gluon four-point vertex 
traced over the adjoint color indices. In this case 

^"'TaMjb^'J.ns) = -g'^mlcpSijf^'" {p,q,r,s) 

= g^cpS^.r''- , (A.46) 

where cf = (iVc - l)/(2iVc). There is no tree-level term. 
The tensor in the HTL correction term is only defined for 

p + q — r + s — 



f'""{p,q,r,s) 



11 V 

y y 



— + — 

s-y 



[{r~p)-y] [{s+p)-y] 



(A.47) 



This tensor is symmetric in fj, and u and is traceless. It satis- 
fies the Ward identity: 

p^r'"'(p,g,r,s)=r"'(g,r-p,s)-r(g,r,s+p) . (A.48) 



6. HTL Quark Counterterm 

The Feynman rule for the insertion of an HTL quark coun- 
terterm into a quark propagator is 



ir''E(p) 



(A.49) 



where E(p) is the HTL quark self-energy given in lA.39t 
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7. Imaginary-time formalism 



APPENDIX B: SUM- INTEGRALS 



In the imaginary-time formalism, Minkoswski energies iiave 
discrete imaginary values po — i{2TinT) and integrals over 
Minkowski space are replaced by sum-integrals over Euclidean 
vectors (27rnT, p). We will use the notation P — (Po,p) for 
Euclidean momenta. The magnitude of the spatial momen- 
tum will be denoted p = |p| , and should not be confused with 
a Minkowski vector. The inner product of two Euclidean vec- 
tors is P ■ Q = PoQo + P ■ q- The vector that specifies the 
thermal rest frame remains n = (1,0). 

The Feynman rules for Minkowski space given above can 
be easily adapted to Euclidean space. The Euclidean tensor 
in a given Feynman rule is obtained from the correspond- 
ing Minkowski tensor with raised indices by replacing each 
Minkowski energy po by iPo, where Po is the corresponding 
Euclidean energy, and multipying by — i for every index. 
This prescription transforms p — (po,p) into P = (Po,p), 
g'"' into -5^", an d p-q i nto -P Q. The ef fect on the HTL 
tensors defined in IIA.2^ . ljA.43^ . and l|A.47|l is equivalent to 
substituting p-n — P-7V where A'^ — (— i,0), p-y — > —P Y 
where Y = (— i,y), and Y^. For example, the Eu- 

clidean tensor corresponding to ljA.2|l is 



In the imaginary-time formalism for thermal field theory, 
the 4- momentum P = (Po, p)is Euclidean with P^ = Pq -|-p^. 
The Euclidean energy po has discrete values: Po = 2miT for 
bosons and Po = (2n + 1)tiT for fermions, where n is an inte- 
ger. Loop diagrams involve sums over Po and integrals over p. 
With dimensional regularization, the integral is generalized to 
d = 3 — 2e spatial dimensions. We define the dimensionally 
regularized sum-integral by 



i 

j{p} 



47r 



~l 2 

47r 



Po = 2n7rT 
Po = {2n + l)7rT 



d'-^^p 

(27r)3-2^ ' 



(2^) 



(B.l) 



(B.2) 



where 3 — 2e is the dimension of space and fi is an arbitrary 
momentum scale. The factor {e'^ /inY is introduced so that, 
after minimal subtraction of the poles in e due to ultraviolet 
divergences, /i coincides with the renormalization scale of the 
MS renormalization scheme. 



T^''(P,-P) = (Y^Y 



. P-N 
P-Y 



(A.50) 



1. One- loop sum-integrals 



The average is taken over the directions of the unit vector y. 

Alternatively, one can calculate a diagram by using the 
Feynman rules for Minkowski momenta, reducing the expres- 
sions for diagrams to scalars, and then make the appropri- 
ate substitutions, such as p^ —P^, P ■ <1 ^ —P ■ Q, and 
n - p ^ in - P. For example, the propagator functions ||A.21|I 
and JA.22II become 



At(P) = 
Al(P) = 



P2+Ut{P) 
1 

p2 + nL(P) ■ 



(A.51) 
(A.52) 



The expressions for the HTL self-energy functions I1t{P) and 
II_l(P) are given by l|A.13fl and l|A.14fl with rip replaced by 



: p /P and T™(p, -p) replaced by 



.(.) 



dc (1 



iPo 



iPo - pc 



(A.53) 



Note that this function differs by a sign from the 00 compo- 
nent of the Euclidean tensor corresponding to 1A.21 : 



T''\P,~P) = ~T''°{p,-p) 



(A.54) 



A more convenient form for calculating sum-integrals that 
involve the function Tp is 



Tp 



p2 
^0 



Ptf + p'^c^ 



(A.55) 



where the angular brackets represent an average over c defined 
by 

{f{c)l = w{e) f dc(l-c2)-7(c) (A.56) 
Jo 



and iu(e) is given in 1A.16II 



The simple one-loop sum-integrals required in our calcula- 
tions can be derived from the formulas 



i 



j> 

(P2)" 



/ fj. 2r(f +m-e)r(n-|-m-he 
Vi^y r(n)r(2 - 2e) 

X r(l - e)C(2n - 2m - 3 + 2e)e'^ 

r,A+2m-2n In \l + 2m-2n 



%P} ip' 



(2- 



XT" 

'2n — 2m — d 



\2n) 



1)] 



(P2)" 



(B.3) 
(B.4) 



The specific bosonic one-loop sum-integrals needed are 



i 



P2 12 UttT; 



l+|2 + 2|lzi^,. 



' 4 +\(-i) ^ C(-i) 



p 

(P2)2 
1 

(p2)2 



ly2 

8 



(47r)2 V47rr/ 



i-f27 + f ^-471 I e 



1 

p2p2 



(47r)2 V47rry 



- + 2-f + 2 



4- I 4 + 47 -I- — - 471 ) e 



. (B.5) 
(B.6) 

(B.7) 

(B.8) 



The specific fermionic one-loop sum-integrals needed are 



(B.9) 
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4 ± ^ ^llfJ^ 



x|l+ ( 2-21og2 + 2|^ le 
+ ^4+^-41og2-21og''2 

+4(1 -1082)04^ + 20^^.^' 



C(-i) C(-i) 



2')2 
2 



+ 27 + 41og2 



2\2 



(47r)2 V47rTy 
16 V47rTy 



l+l^-21og2 + 2i:(^,e 



{P} 



i 

j{p} 



(p2)3 



(p2)3 



(47r)2 V47rTy 
3 



i + 27- ^+41og2 
e 3 



5r^ / M y 

64 K^ttT) 



l+|ii-21og2 + 2^,. 



(47r)2 V47rTy 
5 



8 



i+27-i^+41og2 
e 15 



(B.IO) 
,(B.ll) 

,(B.12) 

(B.13) 

(B.14) 

(B.15) 



1 



1 



-2 + 27 + 41og2 



P2 (47r)2 V47rTy 

+ ( 4 + 81og2 + 41og''2 + 47(l + 21og2) 



+ — -471 )e 



(B.16) 



The errors are all of one order higher in e than the small- 
est term shown. The number 71 is the first Stieltjes gamma 
constant defined by the equation 

C(l + z) = i + 7 - 71^ + C»(^'') • (B.17) 

z 



2. One-loop HTL sum-integrals 

We also need some more difficult one-loop sum-integrals 
that involve the HTL function defined in (IA.33II . 
The specific bosonic sum-integrals needed are 



"1 



x(-l) 



+ 27 + 21og2 



(B.18) 



2p2 



^{j^Y 

(47r)2 V47rry 



2 log 2 



§2(1+27) +21og22 + ^ 



i 



(P2)2 



Tp = 



.2. ^ 




/ 2 


i + 27 + 1 



(B.19) 
.(B.20) 



The specific fermionic sum-integrals needed are 
1 



(p2)2 



Tp = 



(47r)2 V47rry 



- + 27 + l + 41og2 



{P^ 



p2p2 



Tp 



2 / M 
(47r)2 V47rry 



i 



log2(i + 27) + 51og22 + ^ 



(B.21) 



(B.22) 



p2 p2 
{P} ^ ^0 



Tp = 



^(j^Y 



l + 2(7 + 21og2)i + ^ 



. ^2iTpf 
{P} P ^0 



+4 log^ 2 + 87 log 2 -471 



(J^Y 



(47r)2 \4ttTJ 
log 2 (i + 27) -f-51og^2 

1 ( M Y' 



P2\(Q.y)2/^ (47r)2 V47rry 



- - 1 +27 + 41og2 



(B.23) 



(B.24) 



(B.25) 



The errors are all of order e 

It is straightforward to calculate the sum-integrals ljB.21fl - 
IIB.2411 using the representation lA.SSt of the function Tp. For 
example, the sum-integral (IB.18t can be written 



p2 



(B.26) 



I p f J p f \Po ^ P^'^^ 

where the angular brackets denote an average over c as defined 



tTp 



2 2 
p c 



P§ + p2c2 



(B.27) 



The first term in the square brackets vanishes with dimen- 
sional regularization, while after rescaling the momentum by 
p ^ p/c, the second term reads 



Evaluating the average over c, using the expression l|B.8|l for 
the sum-integral, and expanding in powers of e, we obtain 
the result 1B.18L Following the same strategy, all the sum- 
integrals l)B.21^ - l|B.24^ can be reduced to linear combinations 
of simple sum-integrals with coefficients that are averages over 
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c. The only difficult integral is the double average over c that 
arises from IIB.2411: 



= 2 log 2 + 2 (log^ 2 - 2 log 2) e 



(B.29) 



3. Simple two-loop sum-integrals 



The simple two-loop sum-integrals that are needed are 
1 

{PQ} 



= 0, 



i 



„2 



i+4-21og2-f4^i-i^ 

e C(-i) 



{PQ} P2Q2^4 (4^)2 V47rT 



{PQ} 



J{PQ} 









-!- - 7.002 




e 



^ ^ ^ V 

7r)2 V47rT/ 



(4- 



i 



1-67 + 



x'(-i) 



C(-i) 



tpQ} q^P^Q^R^ (4- 



7r)2 V47rTy' 

(I) 



-h 9.5667 



-PQ} 



(B.30) 



(B.31) 



(B.32) 



(B.33) 



(B.34) 



(B.35) 



- 8.1420 



(B.36) 



where R = ~{P + Q) and r = |p-|-qj. The corrections are all 
of order e. To motivate the integration formula we will use 
to evaluate the two-loop sum-integrals, we first present the 
analogous integration formula for one-loop sum-integrals. In 
a one-loop sum-integral, the sum over Po can be replaced by 
a contour integral in po ~ ~iPo'- 



i 



F{P) = hm / ^ / [F{-ipo,p) - F{0,p) 



(B.37) 



where n{po) — l/{e^^° — 1) is the Bose-Einstein thermal dis- 
tribution and the contour runs from —00 to -l-oo above the 
real axis and from -l-oo to —00 below the real axis. This for- 
mula can be expressed in a more convenient form by collapsing 
the contour onto the real axis and separating out those terms 
with the exponential convergence factor n(|po|). The remain- 
ing terms run along contours from — 00 ± ze to and have 
the convergence factor 6^^°. This allows the contours to be 
deformed so that they run from to ±ioo along the imagi- 
nary Po axis, which corresponds to real values of Po = ~ipo- 
Assuming that F{—ipo, p) is a real function of po, i.e. that it 
satisfies P(— ipo,p) = F(— ipo,p)*, the resulting formula for 
the sum-integral is 



^P(P) = j^F{P) 



+ I e(po)n(|pol)2ImP(-ipo + e,p) , (B.38) 
J p 

where e(po) is the sign of po. The first integral on the right side 
is over the (d + l)-dimensional Euclidean vector P = (Po, p) 
and the second is over the (d + l)-dimensional Minkowskian 
vector p = (po,p). 

The two-loop sum-integrals can be evaluated by using a 
generalization of the one-loop formula IB. 3811 : 



J{PQ} 



PQ 



F{P)G{Q)H{R) = / F{P)G{Q)H{R)- / e(po)np(|po|) 2ImP(-2po -h e, p) Re / G{Q)H{R) 



Po = -ipo+e 



e(po)np(|po|)2ImG(-ipo+£,p)Re / H{Q)F{R) 



Po = -ipo-l-s 



+ / e(po)ns(|po|)2Imi/(-ipa + e,p)Re / F{Q)G{R) 
Jp Jq 



Po^-ipo+e 

+ I e(po)np(|pol)2ImP(-ipo -f e,p) / e(go)np(lgol) 2 ImG(-igo + £, q) Re-H"(P) 



e(po)np(jpoj) 2ImG(-ipo +£,p) / e(go)ns(|go|) 2 Im_f/'(-i(jo -f- e, q) ReP(P) 

J q 

e{po)nB{\po\) 2ImH{~ipo + e,p) / e(go)np(|go|) 2ImF(-igo + e,q) ReG(P) 
p •'q 



Ro^i(po+'io)+^ 



(B.39) 
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This formula can be derived in 3 steps. First, express the 
sum over Pq as the sum of two contour integrals over po, 
one that encloses the real axis Impo = and another that 
encloses the line Impo = — Imqo- Second, express the the 
sum over qo as a contour integral that encloses the real-50 
axis. The resulting terms can be combined into the expression 
IB. 3911 . The integrals of the imaginary parts that enter into 
our calculation can be reduced to 



- 1 =!?5i:/(«^+^^p). 



f{-ipo + £,p) 



P 2 
P ± 



(B.40) 



e(po)n(|po|)2ImTp 



f{-ipo + £,p) 



E (c"'+'V(±ip + e,p/c)l . (B.41) 



The latter equation is obtained by inserting the expression 
(IX35ll for Tp, using IB. 4011 . and then making the change of 
variable p — ^ p/c to put the thermal integral into a standard 
form. 

As a simple illustration, we apply the formula IB. 3911 to the 



sum-integral l|B.31^ . The nonvanishing terms are 
1 



{PQ} 



P2Q2J.2 



+ 



n^dpol) 27r(5(po - 



"■F(bol) 27r(5(po -P 



Q2r2 



X / nF{\qo\)2n5(qo -9^)^ 
'q 



(B.42) 



The delta functions can be used to evaluate the integrals 
over po and 50- The integral over Q is given in 10.1121 up to 
corrections of order e. This reduces the sum-integral to 



npip) 

1 

P 



-!- + 4 - 2 log 2 



nF{p)nF{q) 1 
pq 



(B.43) 



The momentum integrals are evaluated in IIC.3t and 1C.41 . 
Keeping all terms that contribute through order e", we get 



the result IB. 3111 . The sum-integral IIB.32t can be evaluated 
in the same way: 



4 

J{PQ} 



(47r)2 



- - 2 log 2 



+ 



nF{p)nF{q) q^ 
pq 



n-F(p) -2e 
P 



(B.44) 



The sum-integral <B.34t can be reduced to a linear combina- 



tion of l|B.31fl and l|B.32fl by expressing the numerator in the 
form P Q = PoQo + (r^ -p^~ q^)/2 and noting that the PoQo 
term vanishes upon summing over Po or Qo- 

The sum-integral l|B.33^ is a little more difficult. After 
applying the formula (IB.39t and using the delta functions to 
integrate over po, go, and ro, it can be reduced to 

ns(p) 



i 



P'^Q^r'^R^ 



npip) 



J.2 „2 



+ 



nF{p)nF{q) p^ - p^ - q^ 
pq r2 A(p,Q,r) 



nF{p)nB{q) 
pq 



2 2 
p r 

q^ q^ 



A{p,q,r) 



(B.45) 



where A(p, q, r) is the triangle function that is negative when 
p, q, and r are the lengths of 3 sides of a triangle: 

A(p, g, r) = p" + + r* — 2{p'^q + qr^ + r^p') . (B.46) 



After using IIC.118t - JU.1201 to integrate over Q, the first term 
on the right side of IIB.45t is evaluated using IIU.3t . The 2- 
loop thermal integrals on the right side of IIB.45t are given 
in IIC.8t - IIC.lH . Adding together all the terms, we get the 
final result <B.33t . The sum- integrals IIB.35t and lB.36t are 
evaluated in a similar manner. 



4. Two-loop HTL sum-integrals 



We also need some more difficult two-loop sum-integrals 
that involve the functions Tp defined in IIA.3311 



1 



-PQ} 



PQ} 



p2Q2j.4 



■Tp, 



Tp = 



J{PQ} 



PQ 

P'^Q^r' 



PQ} 



2 2 
r — p 

P'^q'^QlR'^ 



To 



J_ ( ^ V 

(47r)2 V47rry 

(47r)2 V47rr/ 



1 

"48 



576 J 



^ + ('2-f 121og2-h4^^-^ I -!- + 136.362 



^ + f ^ -H521og2-H4 ^,,^ }} ] i-f 446.438 



■Tp = 



rj-i2 



(47r)2 V47rr/ 



C(-i) 



^ + ("4 log 2 + 4 ^^,^ }} ] i + 69.174 



C(-i) 



(47r)2 \AtvTJ 8 



^ + f 2 + 27-1- ^log2-f2 '^^f }} ] i +46.8757 



C(-i) 



(B.47) 
(B.48) 
(B.49) 
(B.50) 
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The errors are all of order e. To calculate the sum- 
integral ljB.47^ . we begin by using the representation ijA.i 
of the function Tr: 



{PQ} 



P2Q2J.2 



4 



{PQ} 



{PQ} 



(B.51) 



The first sum-integral on the right hand side is given 
by IB. Sit . To evaluate the second sum-integral, we apply 
the sum-integral formula (IB.39t : 



{PQ} P^QHR'o + r^<^- 



-2Re 



+ 



p 

-3-|-2e 



ub (p) f 1 



Pq — — Zp + £ 



P^(-ip,p/c) 



-2c 



pq 

■i+2e 



pq 



A(p + ie, q, rc) 



nF{p)nB{q) 



Re 



2 2 

-q 



pq 



A{p + ie,q,rc) 



(B.52) 



where rc = jp + q/c|. In the terms on the right side with a 
single thermal integral, the appropriate averages over c of the 
integrals over Q are given in ljC.116^ and l|C.124^ . 

The subsequent integrals over p are special cases of llC.311 
and ltCl4|l : 



ns(p)p 



nF{p)p 



= 2= 



(1)-4.(|)2. C(-l + 4e) 
(1)_2.(|)-. C(-l) 



x(eV)^(47rr)-«^, (B.53) 



= [1 



-1+4e] 



nB{p)p' 



(B.54) 



This yields 



,2 / !^Re 



Rf, + r'^c^ 



Po^-ip + s 



+ 



nsip) 



'Q 



P->(-ip,p/c) , 



(47r)2 \4:TtTJ 



121og2-4^:f— f - + 70.122 



(B.55) 



For the two terms in IB. 5111 with a double thermal integral, 
the averages weighted by c are given in iHTTt and HCITI . 
Adding them to llB.55l l. the final result is 



{PQ} 



(47r)2 V47rr/ 



(5) 



4 - fe- 121og2-4 ^,f ]} ] - -1- 51.9307 



c(-i) y ^ 



Inserting this into IB. Sit , we obtain the final result IB. 4711 . 
The sum-integral ljB.48^ is evaluated in a similar way to 
i.47ll . Using the representation HA.SSt for Tr, we get 



PQ} 



p2Q2^4 



Tr = 



4 

^{PQ} 



{PQ} 
2 



p2Q2^i 



(B.57) 



The first sum-integral on the right hand side is given 
by IIB.32II . To evaluate the second sum-integral, we apply 
the sum-integral formula (IB. 3911 : 



J{- 



{PQ] P^Q^r^Rf, + r^c^ 

TIF (p) 



-Re 



+c 



+ 



' p 

-l + 2e 



P 



2 I 2 

p +q 



Q^r'^{R^ + r^c-2 



PQ — — ip+!: 



nsip) -2 

P 

P JQ 



Q^R^ 



nF{p)nF{q) (f_ r'^c^ 



P^(-ip,p/c) 

P'^q' 



pq 



pq 

i+2e f nF{p)nB{q) p'^ +r, 



A{p + ie, q, rc) 



pq 



xRe 



222 

rc-p -q 

A{p + ie,q,rc] 



(B.58) 



In the terms on the right side with a single thermal inte- 
gral, the weighted averages over c of the integrals over Q are 
given in HC. 122(1 . HC.127^ . and HC.128^ : After using ljB.54fl to 
evaluate the thermal integral, we obtain 
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For the two terms in IIB.58t with a double thermal integral, 
the averages weighted by (? are given in EHJ, (EJU, and 
IIC.24L Adding them to 1B.59II . the final result is 
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To evaluate (IB. 4911 . we use the expression (TTSSll for Tr 

and the identity P Q = {R^ - - Q^)/2 to write it in the 
form 

PQ 



J{PQ} ^ ^ ^ J{PQ} 
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The sum-integrals in the first 3 terms on the right side of 
linjTl are given in lIHTTol . jnUHJ, fOTl . and llOH . The 

last sum-integral before the average weighted by c is given in 
(IB. 5211 . The average weighted by (? is given in IIB.56t . The 
average weighted by c* can be computed in the same way. 
In the integrand of the single thermal integral, the weighted 
averages over c of the integrals over Q axe given in lIC. 11711 
and ljC.126|l : After using ljB.54^ to evaluate the thermal inte- 
gral, we obtain 
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For the two terms with a double thermal integral, the averages 
weighted by c are given in JHTSI and lIH^ . Adding them 
to 1B.62II. we obtain 
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We finally need to evaluate IB. 5011 . Applying lB.39t gives 
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In the terms on the right side, with a single thermal fac- 
tor, the weighted average is given in Eq. HC.129|I . After usmg 
Eq. IB. 5411 to evaluate the thermal integral, we obtain 
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Q^r'2{Rl + r^c'2 
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The terms with two thermal factors are given in 
Eqs. jn^Oj, lltnst and JtHel . Adding them to lIH^GSl . we 
finally obtain IB. 5011 . 



APPENDIX C: INTEGRALS 

Dimensional regularization can be used to regularize both 
the ultraviolet divergences and infrared divergences in 3- 
dimensional integrals over momenta. The spacial dimension is 
generalized to d = 3 — 2e dimensions. Integrals are evaluated 
at a value of d for which they converge and then analytically 
continued to d = 3. We use the integration measure 
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1. 3-dimensional integrals 

We require one integral that does not involve the Bose- 
Einstein distribution function. The momentum scale in these 
integrals is set by the mass m = mo- The one- loop integral 
is 



The error is one order higher in e than the smallest term 
shown. 



2. Thermal integrals 

The thermal integrals involve the Fermi-Dirac distribution 
np (p) = l/(e^f + 1). The one-loop integrals can all be ob- 
tained from the general formula 
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The simple two-loop thermal integrals are 
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We also need some more complicated 2-loop thermal integrals 



that involve the triangle function defined in Eq. ljB.4i 
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The most difficult thermal integrals to evaluate involve 
both the triangle function and the HTL average defined in 
lU36ll . There are 2 sets of these integrals. The first set is 
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The second set of these integrals involve the variable rc = 
lP + q/c|: 
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The simplest way to evaluate integrals like ljC.4|l - l|C.7fl whose 
integrands factor into separate functions of p, q, and r is to 
Fourier transform to coordinate space where they reduce to 
an integral over a single coordinate R: 



fip)9iq)h{r) 



The Fourier transform is 
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f{R) = I e'P-^/(p) , (C.28) 



and the dimensionally regularized coordinate integral is 
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The Fourier transforms we need are 
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If a is an even integer, the Fourier transform 1C.31II is partic- 
ularly simple in the limit d ^ 3: 
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where x = ttRT 

We can use these simple expressions only if the integral over 
the coordinate R in (IC.27II converges for d = 3. Otherwise, 
we must first make subtractions on the integrand to make the 
integral convergent. 

The integrals iCH - iCH can be evaluated directly by 
applying the Fourier transform formula 10.2711 in the limit 
0. 

The integrals IIO.8t - lO.10ll can be evaluated by first av- 
eraging over angles. The triangle function can be expressed 
as 

A(p,g,r) = -4p^g^(l-cos^6i) , (0.34) 
where 9 is the angle between p and q. For example, the angle 
average for 10.811 is 
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After integrating over x and inserting the result into 
the integral reduces to 



np{p) np{q) 
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The integrals over p and q factor into separate integrals that 
can be evaluated using ifO^ . After averaging over angles, the 
integrals l|0.9^ and 10. 10^ reduce to 
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The integral ifoTT^ can be evaluated by using the identity 



2 I 2 

p +q 

A(p, q, r) 



1 



2e 1 



8e p2g2 



(0.39) 



The identity can be proved by expressing the angular averages 
in terms of integrals over the cosine of the angle between p 
and q as in ljO.35^ . and then integrating by parts. Inserting 
the identity 10.391 into lO.llll . the integral reduces to 
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The integral in the first term on the right is given in l|0.5^ . 
while the second term can be evaluated using 10. 3t . 

The integral 10.1711 can be evaluated directly in three di- 
mensions by first averaging over c and x, and then integrate 
the resulting functions numerically over p and q. 

To evaluate the weighted averages over c of the thermal 
integrals in Egs. 10.181 - 1O.20L we first isolate the divergent 
parts, which come from the region p — q — > 0. We write the 
product of thermal functions in the form 
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where s = {p + q)/2. In the difference term, the HTL average 
over c and the angular average over a; = p-q can be calculated 
in three dimensions; 
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The remaining 2-dimensional integral over p and g can then 
be evaluated numerically: 
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The integrals involving the n% (s) term in (IC.4H are diver- 
gent, so the HTL average over c and the angular average over 
a; = p • q must be calculated in 3 — 2e dimensions. The first 
step in the calculation of the (s) term is to change variables 
fi-om p and qto s — {p + q)/2, P — ipq/ (p + g)^, and x = p-q: 
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where s± = s[l ± y/T^^ and r = s[4 - 2/3(1 - x)]'^/^. The 2 
terms inside the average over x come from the regions p > q 
andp < g, respectively. The integral over s is easily evaluated: 
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It remains only to evaluate the averages over c and x and the 
integral over l3. 



The first step in the calculation of the n%{s) term of JC.18II 
is to decompose the integrand into 2 terms: 
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(C.43) The weighted averages over c gives a hypergeometric function: 
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In the -l-g case of IIC.5211 . the ie prescription is unnecessary. 
The argument of the hypergeometric function can be written 
1 — Py, where y = {1 — x)/2. After using a transformation 
formula to change the argument to Py, we can evaluate the 
angular average over x to obtain hypergeometric functions 
with argument /3. For example, the average over x of llU.5211 
is 
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where (a)b is Pochhammer's symbol which is defined in 
lie. 14711 . Integrating over /3, we obtain hypergeometric func- 
tions with argument 1: 




Expanding in powers of e, we obtain 
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In the — g case of <C.52ll . the argument of the hypergeo- 
metric functions can be written (1 — Py)/{1 — P ikie), where 
y = (1 — a::)/2 and the prescriptions +ie and —ie correspond to 
the regions p > q and p < q, respectively. These regions corre- 
spond to the two terms inside the average over x in IIC.4811 . In 
order to obtain an analytic result in terms of hypergeometric 
functions, it is necessary to integrate over P before averaging 
over X. The integrals over P can be evaluated by first us- 
ing a transformation formula to change the argument of the 
hypergeometric function to —P{1 — y)/(l — P) and then us- 
ing the integration formula IC. 15411 to obtain hypergeometric 
functions with arguments y or 1 — y: 
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In the +g case of IIC.6111 . the ie prescription is unnecessary. 
The argument of the hypergeometric function can be written 
1 — Py, where y = {1 — x)/2. After using a transformation 
formula to change the argument to /3j/, we can evaluate the 
angular average over x to obtain hypergeometric functions 
with argument (3. For example, the average over x of ljC.61|l 



After averaging over x, we obtain hypergeometric functions 
with argument 1: 
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Expanding in powers of e and then taking the real parts, we 
obtain 
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To evaluate the subtraction in the integral (IC.46L we use 
the identity q'^ = (r^ + — — 2p • q)/2. The integral with 

— in the numerator is purely imaginary. Thus the real 
part of the integral can be expressed as 
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The first term in Eq. I)C.59|I is decomposed into 2 terms: 
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The weighted averages over c give hypergeometric functions: 
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where {a)t is Pochhammer's symbol which is defined in 
l|C.147fl . Integrating over (3, we obtain hypergeometric func- 
tions with argument 1: 
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Expanding in powers of e, we obtain 
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In the —q case of IIU.6111 . the argument of the hypergeo- 
metric functions can be written (1 — Py)/{1 — /3 ± ie), where 
y = {l — x)/2 and the prescriptions -|-ie and — ie correspond to 
the regions p > q and p < q, respectively. These regions corre- 
spond to the two terms inside the average over x in JC.48II . In 
order to obtain an analytic result in terms of hypergeometric 
functions, it is necessary to integrate over (3 before averaging 
over X. The integrals over /3 can be evaluated by first us- 
ing a transformation formula to change the argument of the 
hypergeometric function to — /3(1 — y)/{l ~ 13) and then us- 
ing the integration formula IC. 15411 to obtain hypergeometric 
functions with arguments y or 1 — y: 
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After averaging over x, we obtain hypergeometric functions 
with argument 1: 
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Expanding in powers of e and then taking the real parts, we 
obtain 
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Inserting the sum of the integrals <C.66ll and (IU.69t into the 
thermal integral 1C.48L we obtain 
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It remains only to evaluate the integral in Eq. (IC.59II with 
p ■ q in the numerator. We begin by using the identity 
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In the first term on the right side, the average over c is a 
simple multiplicative factor: {c^)c = 1/(3 — 2e). The average 
over X gives hypergeometric functions of argument /3: 
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3- 2e 
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The integral over /3 gives hypergeometric functions of argu- 
ment 1: 
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Expanding in powers of e, we obtain 
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In the second term of IIC.71L the average over c is given by 
'.621 . In the +q term, the average over x = p ■ q is 



5 - 2e 
4e 



F 



2-e,l,f 
3-2e, 1+e 



-F 



F 



1-6, l,f 
3-2e, 1 + e 

1 - 2e, f - e 
3-3e 



5 (l).(l)-2.(3)-2.(|)-. 

+ 4e (l)-.(3)-3. ^ 



/3 -i^F 

/ 1 - e I 3 - 3e 



(C.75) 



Integrating over f3, we obtain hypergeometric functions of ar- 
gument 1: 
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Expanding in powers of e, we obtain 
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In the —q term in the integral of the second term of ljC.71^ . 
we integrate over /3 before averaging over x. The integral over 
/3 can be expressed in terms of hypergeometric functions of 
type 2-Fi: 
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The phase in the last term is e""^'' for the /(s+,s_,r) term 
of lU.48t , which comes from the p > q region of the integral, 
and e"^'' for the /(s_ , s+, r) term, which comes from the p < q 
region. The average over a; = p • q can be expressed in terms 
of hypergeometric functions of type 3F2 evaluated at 1: 
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The expansion of the real part of the integral in powers of 6 
is 
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Inserting (IC.74II . HO. 7711 . and IIC.8011 into the thermal integral 
of 1IU.7IL we obtain 
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Inserting this along with lU.70t into IIC.5911 . we obtain 
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Adding this integral to the subtracted integral in l)C.46^ . we 
obtain the final result in IC.IQII . The subtracted integral ap- 
pearing in lie. 471 vanishes due to antisymmetry of the inte- 
grand. Thus the final result ljC.20fl is given by HC.47fl . 

The integrals iClTl and lIH^ can be computed directly 
in three dimensions, as described above. The integrals 1C.23II - 
IIU.26t are divergent and require subtractions to remove the 
divergences. We first isolate the divergent part which come 
from the region q 0. We need one subtraction: 



nsiq) 



T 1 \ T 1 
nB(.)--+2)+--2. (C.83) 



In the integral lC.24t , it is convenient to first use the identity 
— p^ + 2p ■ q/c + q^ /c^ to expand it into 3 integrals, two 

of which are ljC.21fl and ljC.23^ . In the third integral, the 

subtraction 1U.83I I is needed to remove the divergences. 
For the convergent terms, the HTL average over c and the 

angular average over a; = p ■ q can be calculated in three 

dimensions: 
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The remaining 2-dimensional integral over p and g can then 
be evaluated numerically: 
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The integrals involving the terms subtracted from n{q) 
in lie. 8311 are divergent, so the HTL average over c and the 
angular average over a; = p • q must be calculated in 3 — 2e 
dimensions. The first step in the calculation of the subtracted 
terms is to replace the average over c of the integral over q by 
an average over c and x: 
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The integral over p can now be evaluated easily using either 
.5411 or 



nAp)p-'-'' = ^^(l-2*^)C(l-4.)(eV)^ri-*^ • 

(C.92) 

It remains only to calculate the averages over c and x. The 
averages over x give 2F1 hypergeometric functions with argu- 
ment [(1 =F c)/2 - ie]"^: 
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Using a transformation formula, the arguments can be 
changed to (l=pc)/2 — ie. If the expressions (IC.93II and 10.9411 
are averaged over c with a weight that is an even function of 
c, the + and — terms combine to give 3F2 hypergeometric 
functions with argument 1. For example, 
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Upon expanding the hypergeometric functions in powers of e 
and taking the real parts, we obtain 
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If the expressions lj(T93|l and lIcTM)! are averaged over c 
with a weight that is an odd function of c, they reduce to 
integrals of 2F1 hypergeometric functions with argument y. 
For example. 
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The resulting expansions for the real parts of the averages 
over c and x are 
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(C.102) 
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Multiplying each of these expansions by the appropriate fac- 
tors from the integral over q in l|C.91^ and the integral over p 
in IIU.92t or IIB.541 1. we obtain 
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Adding Eq. 1C.104II to the subtracted integral 
obtain the final result in Eq. 1C.23II . Combining 
IIC.105II and ilU.10611 . we obtain 
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The integral lIHlH is obtained from iClTl . 

and He. 1071 1. Finally consider lIH^ and (IHlel . In order 

to evaluate them we need two subtractions for each integral 
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The subtractions can be evaluated directly in three dimen- 
sions and the results are given in Eqs. (IC.89t - l|C.90t The in- 
tegrals He. 2511 and <U.26t are then given by the by the sum 
of the difference terms (IU.89II and <C.90t and the subtraction 
terms HC.lOSII - IIC.lllll . 



3. 4-dimensional integrals 



In the sum-integral formula (IB.39II . the second term on the 
right side involves an integral over 4-dimensional Euclidean 
momenta. The integrands are functions of the integration 
variable Q and R = —{P + Q). The simplest integrals to 
evaluate are those whose integrands are independent of Pq: 
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Another simple integral that is needed depends only on P 
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where {a)b is Pochhammer's symbol which is defined in 
IC. 14711 . We need the following weighted averages over c of 
this function evaluated at P = (— ip, p/c): 
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The remaining integrals are functions of Pq that must be 
analytically continued to the point Pq = —ip + e. Several of 
these integrals are straightforward to evaluate: 
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We also need a weighted average over c of the integral in 
1C.118II evaluated at P = (— jp, p/c). The integral itself is 
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The most difficult 4-dimensional integrals to evaluate in- 
volve an HTL average of an integral with denominator + 
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The analytic continuation to Pq — —ip + e is implied in these 
integrals and in all the 4-dimensional integrals in the remain- 
der of this subsection. 



We proceed to describe the evaluation of the integrals 
ijC. 124(1 and HC.126fl . The integral over Qo can be evalu- 
ated by introducing a Feynman parameter to combine and 
Ro + r^c? into a single denominator: 
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where we have carried out the analytic continuation to Po = 
— ip 4- £. Integrating over r and then over the Feynman pa- 
rameter, we get a 2 Pi hypergeometric function with argument 



1-c^ 



Q2(p2 + r2c2) (47r^2 

,..(1)-2.(1) 



2 (e> ) P 



2e, 1 - e 



(2)_3. (^-^^)-^M%-3. 



1 ^ 

1 — c 



(C.131) 



The subsequent weighted averages over c give 3P2 hypergeo- 
metric functions with argument 1: 
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After expanding in powers of e, the real part is 10.1261 . 

The integral <C. 12711 has a factor of 1/r^ in the integrand. 
After using (lU. 13011 . it is convenient to use a second Feyn- 
man parameter to combine (1 — a; -|- a;c^)r^ with the other 
denominator before integrating over r: 
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After integrating over r and then y, we obtain 2P1 hypergeo- 
metric functions with arguments x{l — c^). The integral over 
x gives a 2 Pi hypergeometric function with argument 1 — c^: 
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After averaging over c, we get a hypergeometric functions with 
argument 1: 



1 



(47r) 



/ 7 2^e -2-2e (l)e 



1 (-i)-.(l)-. 1 _ (-i)_.(l)_2.(2)-2. 



2^ 



/ 1 - 2e i - 2e,-e 



(|)-2.(l)-3. 

lU . (C.136) 



Q2r2 \ i?2 +r2c2 / (47r 



r(e^M ) P 



(|)-2e 



(1)-2.(1)-. 

^ 2 ' 



xF 



(|)-2.(l)-3, 

1 - 2e, i - 2e, -e 
I - 2e, -3e 



1 



(C.137) 



After expanding in powers of e, the real part is IIC.127t . 

To evaluate the integral ljC.128^ . it is convenient to first 
express it as the sum of 3 integrals by expanding the factor 
of 5^ in the numerator a,s = + 2p ■ r + r^: 
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To evaluate the integral with p • r in the numerator, we first 
combine the denominators using Feynman parameters as in 
lU. 13411 . After integrating over r and then y, we obtain 2F1 
hypergeometric functions with arguments x{l — c^). The in- 
tegral over X gives 2F1 hypergeometric functions with argu- 
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After averaging over c, we get a hypergeometric function with 
argument 1: 
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After expanding in powers of e, the real part is 
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Combining this with (lU. 12411 and IIC.126L we obtain the in- 
tegral JcTT28l. 

To evaluate the integral (IC.129II . we first express the nu- 
merator as a sum of two integrals whose averages have been 
calculated: 
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The two hypergeometric functions are now combined into a 
single hypergeometric functions, which yields 
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1 (1)41),.(|)-. 
(|)-2. 



-1 + e 



i^e (l)^-2e 



(C.144) 



Expansion in powers of e, yields Eq. IC. 12911 . 



4. Hypergeometric functions 

The generalized hypergeometric function of type pFq is an 
analytic function of one variable with p + q parameters. In 
our case, the parameters are functions of e, so the list of pa- 
rameters sometimes gets lengthy and the standard notation 
for these functions becomes cumbersome. We therefore intro- 
duce a more concise notation: 



F 



Qi, a2, • . • , ftp 

Pi,. . . ,(3g 



z = pFg(ai,a2, . . . ,Qp;/3i, . . . ,/3,; z) 



(C.145) 

The generalized hypergeometric function has a power series 
representation: 



F 



Qi, a2, ■ ■ ■ ,cep 



{ai)n{a2)n ■ ■ ■ {ap)„ „ 



E[ai )n[a2 )n ■ ■ 



(C.146) 
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where {a)b is Pochhammer's symbol: 

r{a + b) 



{a)t 



r(a) 



(C.147) 



The power series converges for \z\ < 1. For z — 1, it converges 
if Res > 0, where 



(C.148) 



The hypergeometric function of type p+iFg+i has an integral 
representation in terms of the hypergeometric function of type 



\ Pi , • ■ • , p? 



tz 



r(^)r(i/) ^ I ai,a2,. . . ,ap,u 



(C.149) 



If a hypergeometric function has an upper and lower param- 
eter that are equal, both parameters can be deleted: 



^ I ai,a2, . . . ,ap,!/ 



= F 



Qi, a2, . . . , Qp 



(C.150) 

The simplest hypergeometric function is the one of type 
\Fo. It can be expressed in an analytic form: 



lJo(c 



= il-z) 



(C.151) 



The next simplest hypergeometric functions are those of type 
2-Fi. They satisfy transformation formulas that allow an 2^1 
with argument z to be expressed in terms of an 2F1 with 
argument z/{z — 1) or as a sum of two 2^1 's with arguments 
1 — 2; or l/jz or 1/(1 — 2). The hypergeometric functions of 
type 2F1 with argument z = 1 can be evaluated analytically 
in terms of gamma functions: 




r(/3i)r(/3i 



a\ — a2} 



r(/3i - Qi)r(/3i - 02) 



(C.152) 



The hypergeometric function of type zF2 with argument z = 1 
can be expressed as a 3-F2 with argument 2 = 1 and different 
parameters [39l |: 



F 



Ql, a2, Ot-i 
Pi, 132 



1 = 



xF 



r(/3i)r(/32)r(s) 
r(ai + s)r(«2 + s)r(a3 



f3i — a3 , P2 — ct3 , s 
ai + s,a2 + s 



(C.153) 



where s = /3i + /32 — ai — a2 — as • If all the parameters of a 3 
are integers and half-odd-integers, this identity can be used to 
obtain equal numbers of half-odd-integers among the upper 
and lower parameters. If the parameters of a 3^2 reduce to 
integers and half-odd-integers in the limit e ^ , the use of 
this identity simplifies the expansion of the hypergeometric 
functions in powers of e . 



The most important integration formulas involving 2F1 hy- 
pergeometric functions is <C.149t with p = 2 and ? = 1. An- 
other useful integration formula is 




r(Qi+M)r(Q2 + M)r(/?i)r(-M) 
r(Qi)r(a2)r(/3i + ^i) 



xF 



Ql -I- /i, Q2 + /i, J' + M 




This is derived by first inserting the integral representation 
for 2^1 in IC. 14911 with integration variable t' and then eval- 
uating the integral over t to get a 2F1 with argument 1 -I- t'z. 
After using a transformation formula to change the argument 
to —t'z, the remaining integrals over t' are evaluated using 
lU. 14911 to get 3i<2's with arguments —z. 

For the calculation of two-loop thermal integrals involving 
HTL averages, we require the expansion in powers of e for 
hypergeometric functions of type pFp-i with argument 1 and 
parameters that are linear in e. If the power series representa- 
tion lU. 14611 of the hypergeometric function is convergent at 
z = 1 for e = 0, this can be accomplished simply by expanding 
the summand in powers of e and then evaluating the sums. If 
the power series is divergent, we must make subtractions on 
the sum before expanding in powers of e. The convergence 
properties of the power series at z = 1 is determined by the 
variable s defined in IIC.148t . If s > 0, the power series con- 
verges. If s — > in the limit e ^ 0, only one subtraction is 
necessary to make the sum convergent: 



F 



Ql, a2, . . . , Op 

f3i,. . . , f3p-i 



r(/?i).-.r(/3p-i) 
r(ai)r(a2)---r(Qp) 



C(s + i) 



+ 



E 



(iai)n(a2)ra ■ ■ ■ (Qp)n 

(/3i)„---(/3,)„n! 

r(/3i)---r(/3p-i) 
r(ai)r(a2)---r(a. 



(C.155) 



If s — > —1 in the limit e — > 0, two subtractions are necessary 
to make the sum convergent: 



F 



Ltl , Lt2 , ■ . ■ , Ltp 

/3i,...,/3p-i 



r(A)---r(/3p,i) 
r(Qi)r(Q2)---r(ap) 



[C(s + i)+K(s + 2)] + 



(ai)„(a2)n ■ • ■ (ap)n 
(/3i)„---(/3g)„n! 



,ai r(a2)---r(ap) ^ ^ J 



n 



where t is given by 



^-E 

i=l 



{ai - l)(a, -2) 



E 



(ft - l)(ft - 2) 



(C.156) 



(C.157) 
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The expansion of a pFp-i hypergeometric function in pow- 
ers of e is particularly simple if in the limit e ^ all its pa- 
rameters are integers or half-odd-integers, with equal numbers 
of half-odd-integers among the upper and lower parameters. 
If the power series representation for such a hypergeometric 
function is expanded in powers of e, the terms in the sum- 
mand will be rational functions of n, possibly multiplied by 
factors of the polylogarithm function ^p{n + a) or its deriva- 
tives. The terms in the sums can often be simplified by using 
the obvious identity 



evaluated using 



J2^f{n)-f{n + k)] = ^/(i) 



(C.158) 



The sums over n of rational functions of n can be evaluated 
by applying the partial fraction decomposition and then using 
identities such as 

oo 

oc 

The sums of polygamma functions of 71 -I- 1 or n -f- ^ divided 
by n -I- 1 or n -I- i can be evaluated using 



E 



E 

n=0 





log(n + 1) 


n + 1 


n + 1 


tpin + l) 


log(n -1- 1) 


n+i 


n+1 




log(n + 1) 


n+1 


n + 1 




log(n + 1) 


n+i 


n + 1 



'27 -J^-f^ 



(C.161) 



+ — -71 , (C.162) 
= -il" -41og2-^21og'2 



IT 
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71 , (C.163) 



--(7 + 21og2)^ 



-— -71 , (C.164) 

where 71 is Stieltje's first gamma constant defined in <B.17t . 
The sums of polygamma functions of n -f 1 or n -I- i can be 



oc y \ 

EU(-+i)-i°g("+i) + 2(;^) 

n=0 ^ ^ ' / 



= ^ + ^7-ilog(2^) , (C.165) 



E (^(^ + I) - log(" + 1) + 



i7-log2- ilog(2^). (C.166) 



We also need the expansions in e of some integrals of 2F1 
hypergeometric functions of y that have a factor of |1 — 2y|. 
For example, the following 2 integrals are needed to obtain 



dyy-''{l-yy+'\l-2y\F 



l-e,e 
-3e 



i + (^ + ^log2)., (C.167) 



dyy'+^il-yf+^\l-2y\F 



2 + 2e, 1 + e 
2 + 3e 



J+(^ + |log2)e. (C.168) 



These integrals can be evaluated by expressing them in the 
form 



dyy-'\l-yr-^\l-2y\F 



dyy^-\l^yr-'{2y^l)F{ 

v ^1 



+ 2 dyy''-\l-yr-\l-2y)F 



Ql, Q2 

/3i 



(C.169) 



The evaluation of the first integral on the right side gives 
3F2 hypergeometric functions with argument 1. The integrals 
from to ^ can be evaluated by expanding the power series 
representation IIC.146t of the hypergeometric function in pow- 
ers of e. The resulting series can be summed analytically and 
then the integral over y can be evaluated. 
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